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State-to-state  reaction  probabilities  are  found  to  be  highly  final-state  specific  at  state-selected 
threshold  energies  for  the  reactions  O  Hj  —OH  -i-  H  and  H  -(-  Hj  —  Hj  -(-  H.  The  study 
includes  initial  rotational  states  with  quantum  numbers  0-lS,  and  the  specificity  is  especially 
dramatic  for  the  more  highly  rotationally  excited  reactants.  The  analysis  is  based  on  accurate 
quantum  mechanical  reactive  scattering  calculations.  Final-state  specificity  is  shown  in  general 
to  increase  with  the  rotational  quantum  number  of  the  reactant  diatom,  and  the  trends  are 
confirmed  for  both  zero  and  nonzero  values  of  the  total  angular  momentum. 


I.  INTRODUCTION 

Quantum  mechanical  scattering  theory  allows  us  to 
study  many  characteristics  of  state-to-state  dynamical  pro¬ 
cesses  more  systematically  and  at  a  greater  level  of  detail 
than  is  currently  feasible  by  experiment.  A  classic  example  is 
the  detailed  analysis  of  statistical  vs  state-specific  behavior  in 
individual  elements  of  the  scattering  matrix  (S  matrix)  for 
inelastic  scattering  by  Lester  and  Bernstein. '  At  the  level  of 
the  S  matrix,  we  can  observe  transition  probabilities  between 
individual  channels  specified  completely  by  their  full  sets  of 
quantum  numbers — the  ultimate  level  of  resolution  allowed 
by  quantum  mechanics.  The  study  of  transition  probabilities 
at  this  most  detailed  level  is  the  eventual  limit  of  any  increas¬ 
ingly  exact  dynamical  theory  of  chemical  processes  and — in 
addition — it  is  sometimes  directly  relevant  to  certain  mod¬ 
em  experiments  involving  laser-prepared  initial  states  and 
state-specific  product  detection.  In  the  present  paper,  we  re¬ 
port  an  interesting  and  qualitatively  unexpected  behavior 
observed  at  the  5  matrix  level  for  reactive  scattering  involv¬ 
ing  reactants  with  varying  levels  of  internal  excitation,  in¬ 
cluding  highly  rotationally  excited  states  whose  reactions 
are  hard  to  study  experimentally. 

We  have  found  that  state-selected  reactivity  is  charac¬ 
terized  by  a  surprisingly  high  degree  of  final-state  specificity 
at  energies  near  state-selected  reaction  thresholds  for  the  re¬ 
actions  O -h  — OH -t- H  and  H -f  Hj  —  Hj -I- H.  In 
many  cases,  we  find  that  50%  or  more  of  the  state-selected 
reactive  flux  is  into  three  or  fewer  final  states  and  75%  or 
more  is  into  six  or  fewer  final  states,  even  when  up  to  78  final 
states  are  open.  The  trends  have  been  confirmed  for  both 
zero  and  nonzero  values  of  the  total  angular  momentum. 
This  is  an  exciting  result  because  such  a  high  degree  of  final- 
state  specificity  is  unexpected,  particularly  at  high  energy. 

II.  CALCULATIONS 

The  reactive  scattering  calculations  for  O  -I-  H^  were 
carried  out  with  the  realistic  Johnson-Winter-Schatz 


(JWS)  potential  energy  surface,^  and  the  calculations  for 
H  -h  Hi  used  the  highly  accurate  double  many-body  expan¬ 
sion  (DMBE)  surface.^  The  calculations  were  performed 
using  a  linear  algebraic  variational  basis-set  method  based 
on  the  generalized  Newton  variational  principle  for  the  am¬ 
plitude  density.'*  For  both  reactions,  convergence  of  the 
state-to-state  reaction  probabilities  was  estimated  by  com¬ 
paring  the  production  runs  to  convergence  checks  in  which  a 
sufficient  subset  of  basis  set  and  numerical  parameters  was 
varied  to  check  all  sources  of  possible  error  in  the  solution  of 
the  Schrodinger  equation.  The  average  deviation  of  state-to- 
state  reaction  probabilities  grca»*r  than  10  "  *  between  these 
two  runs  was  0.1%-3%,  depending  on  the  system,  energy, 
and  angular  momentum  for  the  results  discussed  in  this  pa¬ 
per.  Detailed  presentations  of  basis  sets  and  numerical  pa¬ 
rameters  for  these  systems  have  been’  *  and  will  be^  present¬ 
ed  elsewhere.  We  simply  mention  that  the  number  of 
coupled  channels  (after  taking  advantage  of  total  angular 
momentum  J,  parity  P,  and  arrangement  symmetry  to  form 
linear-combination  basis  functions  that  are  maximally  de¬ 
coupled*  )  for  the  calculations  in  this  paper  ranges  from  132 
(at  low  energies)  to  196  (at  high  energies)  for  O  Hj  and 
from  162  (for  /  =  0)  to  260  (for  /=  4  and  P—  -t-  1)  for 
H  -I-  Hi  at  all  energies. 

A  state-selected  reaction  probability  P  J!  ( EJ)  for  initial 
state  y,  total  energy  E,  and  total  angular  momentum  J  is 
defined  as  the  sum  of  state-to-state  reaction  probabilities 
P  ^  ( E.J)  over  all  final  states,  where  y  ( y' )  is  a  set  of  diatom 
quantum  numbers  for  reactants  (products) 

Pr  =  ^P%{EJ).  (1) 

The  set  y  consists  of  quantum  numbers  for  arrangement  ( a ) , 
vibration  (u),  and  rotational  angular  momentum  (y),  and 
the  superscript  R  denotes  that  the  final  arrangement  is 
diSereiu  irom  the  initial  arrangement  .  In  both  cases  stud¬ 
ied  here,  there  are  two  identical  final  reactive  arrangements 
(a  =  2  and  3 )  since  Hi  is  homonuclear  and  P  ^  denotes  the 
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reaction  probability  for  either  product  arrangement,  say 
a  =  2.  ( None  of  our  conclusions  would  be  changed  if  we 
summed  over  a  =  2  and  3,  but  it  is  necessary  to  make  a 
consistent  choice  for  tabulation  and  discussion  purposes, 
and  we  chose  to  deal  with  the  individual  a  =  1  -»a’  =  2 
probabilities.)  Thus,  if  5^^  y,(£)  denotes  an  S  matrix  ele¬ 
ment  for  total  energy  E,  total  angular  momentum  J,  parity  P, 
and  initial  and  Anal  orbital  angular  momentum  of  relative 
translation  /  and  / ',  then 


= 


V 


2  min(y,y)  1 


(2) 


where  n  denotes  a  channel,  i.e.,  a  combination  of  values  y„ 
and  l„  of  y  and  l,y(y)  denotes  the  set  of  all  channels  which 
correspond  to  state  y,  and  P  is  a  channel-to-channel  reac¬ 
tion  probability  defined  by 


=  (3) 

Note  that  a  particular  combination  of  y„  and  /„  uniquely 
specifies  the  parity  P,  so  it  is  not  necessary  to  include  a  sum¬ 
mation  over  the  parity  Pon  the  right-hand  side  of  Eq.  (3). 
The  sum  over  n  in  Eq.  (2)  automatically  picks  up  contribu¬ 
tions  from  both  parities. 

For  total  angular  momentum  7  equal  to  zero,  P  is  the 
same  as  since  there  is  only  one  .allowed  .alue  of  /  for 
each  y  and  one  allowed  value  of  / '  for  each  thus  we  can  use 
either  the  state-selected  or  the  channel-selected  language  for 
7  =  0.  In  the  case  of7  f^O,  state-selected  and  channel-select¬ 
ed  results  differ  and  we  will  consider  both;  the  state-selected 
results  are  more  closely  related  to  observables  when  scatter¬ 
ing  angles  are  not  detected,  and  the  channel-selected  results 
are  of  special  interest  because  ( as  stated  in  the  Introduction ) 
they  provide  the  ultimate  level  of  detail  about  reaction  dy¬ 
namics  that  is  allowed  by  quantum  mechanics. 

For  the  purposes  of  this  article,  we  define  a  state-select¬ 
ed  reaction  threshold  as  the  first  maximum  in  the  den¬ 
sity  of  state-selected  density  reaction  probability 
which  is  defined  as  the  energy  derivative  of  P*(£,7),* 


Pr=~P^AEJ),  (4) 

and  similarly  a  channel-selected  reaction  threshold  is  the 
first  maximum  of 

Pn=-~^P^AEJ),  (5) 

dE 

where  P  *  ( EJ') ,  in  analogy  to  P  *  ( £,7) ,  is  a  channel-selected 
reaction  probability  given  by 


p’:,  =  'lPt{EjY  (6) 

n' 

We  will  also  discuss  analogous  quantities  for  final  states. 
We  define  a  final-state-selected  reaction  probability  P  ^  as 

Py='ZP'^iEJ)  (7) 

r 

and  a  density  of  final-state-selected  reaction  probability 
Py(.EJ)  as 

Py=-^P^y(.EA),  (8) 


For  final  states,  £,h,  is  defined  as  the  energy  of  the  first  maxi¬ 
mum  inpy  (Eyf). 

The  present  paper  is  based  on  calculations  for  O  -I-  H, 
with  7  =  0  and  H  -I-  Hj  with  7  =  0  and  7  =  4  over  a  wide 
range  of  energies  (£  =  0.60-1.90  eV  for  O  +  H,  and  0.80- 
1 .60  eV  for  H  -f  Hj ) .  We  examined  several  initial  states  and 
channels  for  each  reaction  and  found  the  thresholds  as  just 
defined  for  the  7  =  0  calculations.  We  determined  the 
threshold  energies  for  H  -I-  Hj  with  7  =  4  by  adding  an  ener¬ 
gy  for  rotation  £7(7  -I-  1 )  =  0.026  eV  to  the7  =  0  energies, 
using  a  value  of  10.6  cm  “  '  =  0.001  3 1  eV  for  B  determined 
from  a  previous  analysis  of  quantized  transition  state  ener¬ 
gies.*  The  7=4  analysis  presented  here  was  carried  out  for 
the  energies  nearest  these  threshold  energies  on  a  preselected 
grid  of  energies  that  span  the  region  of  interest.  The  energy 
used  was  within  0.016  eV  of  the  threshold  energy  calculated 
using  B  in  all  cases.  We  then  examined  the  distributions  of 
final  states  or  channels  in  each  of  the  7  =  0  or  4  threshold 
calculations  in  order  to  examine  the  final-state  specificities 
that  form  the  title  subject  of  this  article. 

The  central  quantities  under  study  here  are  product  dis¬ 
tribution  specificities  for  various  initial  states  or  channels  at 
their  reaction  thresholds.  Thus  we  define  normalized  prod¬ 
uct  distributions  as 

(9) 

Py 

and 

D(v'J.r\vJ,!)=^-  (10) 

P. 

(Since  P^  and  P*„.  include  only  a  particular  a-* a'  pair, 
namely  1  —  2,  the  denominators  are  not  unity. )  We  also  de¬ 
fine  a  partially  summed  distribution 

D{u'\v,j)  =  '^D(v’,f\v,j)  (11) 

/ 

in  order  to  study  the  specificity  of  reaction  into  each  final 
vibrational  level. 

III.  RESULTS  AND  DISCUSSION 

III.  A.  State-selected  specificities  for  O-hH; 

Table  I  summarizes  the  specificity  of  reaction 
(a  =  1  — a'  =  2)  out  of  initial  states  in  the  ground  vibration¬ 
al  manifold  for  the  O  +  Hj  reaction  with  total  angular  mo¬ 
mentum  equal  to  zero.  The  first  four  columns  give  the  value 
of  the  initial  rotational  quantum  number  j,  the  energy  £,hr  of 
the  state-selected  threshold  ( relative  to  the  potential  energy 
of  infinitely  separated  reactants  at  classical  equilibrium ),  the 
relative  translational  energy  £„i,,hr.  and  the  number  of  open 
final  states  N,^  (£)  with  a'  =  2.  More  precisely,  we  define 
the  state-selected  threshold  energy  as  the  first  maxi¬ 
mum  ofp„  (£,7),  and  the  energy  £„i,,hr  is  the  relative  trans¬ 
lational  energy  for  the  given  initial  state  and  total  energy 
£,hr.  To  create  the  table,  the  final  states  with  a'  =  2  were 
ordered  by  decreasing  magnitude  of  the  normalized  state-to- 
state  product  conditional  probabilities  D{v' ,f\v,j).  These 
condi'Trnal  probabilities  were  then  summed,  beginning  with 
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TABLE  I.  Final-state  specificity  of  reaction  out  cf  initial  states  in  the  u  =  0  manifold  for  the  O  -i-  reaction  with  7  =  0. 


j 

F.H, 

(eV) 

(eV) 

25% 

50% 

75% 

mant  (u,/).  D(v'.j\0,j)‘ 

N 

F 

N 

F 

F 

Doir 

0 

0.69 

0.42 

12 

2 

0.17 

3 

0.25 

5 

0.42 

(0.4)  18 

(0.3)  17 

(0,5)  17 

(0.2)  14 

1 

0.69 

0.41 

12 

2 

0.17 

3 

0.25 

5 

0.42 

(0,4)  18 

(0.5)  17 

(0,3)  17 

(0.2)  13 

2 

0.69 

0.38 

12 

2 

0.17 

3 

0.25 

5 

0.42 

(0,5)  18 

(0.4)  17 

(0.6)  16 

(0,3)  14 

3 

0.70 

0,34 

12 

2 

0.17 

4 

0.33 

6 

0.50 

(0.5)  17 

(0,4)  17 

(0.6)  15 

(0.3)  14 

4 

0,71 

0.30 

13 

2 

0.15 

4 

0.31 

6 

0.46 

(0,5)  17 

(0,6)  16 

(0.4)  15 

(0.7)  13 

5 

0.89 

0.41 

23 

2 

0.09 

3 

0.13 

4 

0.17 

(0.9)  23 

(0,10)  23 

(0,8)  16 

(0,11)  15 

6 

0.95 

0.38 

25 

2 

0.08 

3 

0.12 

6 

0.24 

(0.11)  22 

(0,10)  19 

(0,12)  16 

(0,9)  10 

7 

0.96 

0.29 

26 

2 

0.09 

3 

0.12 

7 

0.27 

(0,12)  21 

(('.ID  19 

(0,13)  12 

(0,10)  9 

8 

1.10 

0.33 

30 

2 

0.07 

4 

0.13 

7 

0.23 

(0.14)  21 

(0,13)  18 

(1,8)  11 

(0.15)  9 

9 

1.23 

0.34 

37 

2 

0.05 

3 

0.08 

6 

0.16 

(0,15)  21 

(0,16)  19 

(1.11)  16 

(1.10)  12 

10 

1.38 

0.36 

48 

2 

0.04 

3 

0.06 

6 

0.13 

(0,17)  23 

(1.13)  19 

(0,18)  13 

(0,16)  10 

11 

1.48 

0.32 

52 

2 

0.04 

3 

0.06 

6 

0.12 

(1.14)  21 

(0.18)  17 

(1.15)  16 

(0.19)  14 

12 

1.60 

0.30 

58 

; 

0.02 

3 

0.05 

7 

0.12 

(1.16)  29 

(0.20)  16 

(0,19)  11 

(1.15)  6 

13 

1.72 

0.27 

70 

2 

0.03 

3 

0.04 

8 

0.11 

(1.17)  20 

(1.18)  19 

(0,21)  15 

(1.15)  7 

14 

1.86 

0.25 

78 

1 

0.01 

4 

0.05 

10 

0.13 

(1.19)  31 

(0.22)  10 

(1.17)  6 

(0.23)  6 

‘The  four  final  stales  which  contribute  most  to  F*  and  their  normalized  conditional  probabilities  expressed  as  percentages,  D{v,j\0,j)  x  100%.  listed  in 
order  of  the  size  of  D(i)',f\0,j),  beginning  with  the  greatest  contributor. 


the  largest,  until  the  sum  was  0.25  or  greater.  The  number  of 
contributors  required  to  reach  this  value  of  the  sum  is  listed 
in  the  column  labeled  N(.25%).  Similarly,  the  number  of 
states  contributing  when  the  sum  reached  0.50  and  0.75  are 
shown  in  the  columns  labeled  Ni50% )  and  NO 5%).  Also 
shown  are  columns  headed  F{15%),  F(50%),  and 
/'(75%),  which  denote  what  fraction  of  the  energeti¬ 

cally  accessible  states  contribute  these  percentages  of  the 
normalized  distribution,  i.e., 

F(x%)=^^!^^^.  (12) 

"open 

Finally,  the  four  final  states  which  contribute  most  to  the 
state-selected  reaction  probability  are  listed  at  the  right,  to¬ 
gether  with  their  normalized  conditional  probabilities  ex¬ 
pressed  as  percentages. 

It  is  seen  in  Table  1  that  as  j  is  varied  from  0  to  14,  the 
state-selected  threshold  energy  increases  from  0.69  to  1.86 
eV  and  the  number  of  final  OH  states  open  at  threshold  in¬ 
creases  from  1 2  to  78.  Nevertheless,  for  all  j,  at  most  two  final 
states  are  required  to  account  for  25%  of  the  state-selected 
reactive  flux,  and  for  j  =  12  and  14,  a  single  final  state  ac¬ 
counts  for  over  25%  of  the  flux.  Thus  much  of  the  state- 
selected  reactivity  is  very  final-state  specific,  particularly  at 
high  j.  Perhaps  even  more  dramatic  is  that  in  most  cases, 
three  final  states  account  for  more  than  half  of  the  flux,  and 
about  six  states  account  for  75%  of  the  flux.  The  amount  of 
flux  into  the  single  one  or  two  greatest-contributing  final 
states  becomes  larger  for  high-energy  (i.e.,  high-7)  thresh¬ 
olds,  while  the  remaining  flux  tends  to  be  spread  among 
more  final  states  for  these  thresholds, -as  evidenced  by  the 
larger  values  (7,  8,  and  10)  of  )V(75%)  for  the  three  states 
with  highest  j.  The  fractions  F{x% )  for  all  three  values  of  x 
show  that  the  reactive  flux  becomes  focused  into  a  smaller 
fraction  of  the  open  final  states  as  j  is  increased. 

These  trends  are  demonstrated  in  Figs.  l(a)-l(c),  in 


which  the  value  of  P  ^  at  the  state-selected  threshold  energy 
is  plotted  vs  j  for  initial  states  (0,0),  (0,8),  and  (0,14), 
where  the  first  number  in  parentheses  denotes  v  and  the  sec¬ 
ond  denotes  j.  It  is  seen  that  as  initial  j  is  increased  within  the 
ground  vibrational  manifold,  the  distribution  D{v',f\0,j) 
becomes  peaked  about  a  narrower  range  of  f  values,  and 
vibrationally  nonadiabatic  transitions  become  more  impor¬ 
tant. 

Threshold  reactivity  for  the  v  =  1  manifold  is  summar¬ 
ized  in  Table  II.  The  vibrationally  excited  states  are  even 
more  final-state  specific  than  the  ground  vibrational  states. 
For  four  of  the  12  initial  states  listed  in  Table  II,  a  single  final 
state  accounts  for  over  25%  of  the  reactive  flux.  The  three 
highest-/  states  are  among  these  four,  which  confirms  the 
remarkable  degree  of  final-state  specificity  at  high  /  noted 
above.  This  phenomenon  is  demonstrated  for  the  initial  state 
(1,10)  in  Fig.  2,  in  which  is  plotted  vs  f.  As  for  initial 
state  (0,14)  [Fig.  1(c)],  Z)(i)',/'|1,10)  is  peaked  about  a 
narrow  range  of f  values  for  the  dominant  v',  and  nonadiaba¬ 
tic  transitions  dominate  the  state-selected  reactivity. 

The  larger  values  of  the  product  distribution  specifici¬ 
ties  D{v',f\v,j)  tend  to  be  very  simple  functions  of f,  for  a 
given  t;',  with  a  single  maximum.  One  might  ask  whether  this 
is  3  consequence  of  a  simple  dependence  on  final  internal 
energy.  For  example,  consider  the  four  most  important  final 
states  in  Table  I  for/  =  11, 12,  and  13.  If  the  most  important 
(i>'  =  0,/')  state  were  always  the  one  intermediate  in  energy 
between  the  two  most  important  (v'  =  \,f)  states,  then  we 
would  predict  the  f  values  of  20,  21,  and  22  to  be  the  most 
important  ones  in  the  u’  =  0  manifold  for/  =  1 1 , 1 2,  and  1 3, 
respectively,  but  instead  the  table  shows  that  the  dominant 
(u'  =  0,/)  state  is /  =  18, 20,  and  21  in  these  three  cases.  In 
addition,  some  states  with  very  similar  internal  energies  to 
the  dominant  {v',f)  state  have  very  small  probabilities.  For 
example,  consider  the  /  =  13  case  of  Table  I  again.  In  this 
case,  the  four  final  states  with  the  largest  values  of 
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FIG.  1 .  Nonnalized  product  distribution  D(v',/\p,  j)  and  state-to-state  re¬ 
action  probabilities  vs  final  rotational  quantum  number  /  for  the 
O+Hj  reaction  with  /  =  0  for  (a)  initial  state  ( 0,0)  at  energy  =  0.69 
eV;  (b)  initial  state  (0,8)  at  l.l0eV;and  (c)  initial  state  (0,14)  at  1.86eV. 


i>(o',7’lO,13),  listed  in  Table  I,  have  internal  energies  vary¬ 
ing  between  1.31  and  1.52  eV.  There  are  seven  other  final 
states  with  internal  energies  in  this  mng^  but  the  average 
value  of  ( o',  / 10, 1 3 )  for  these  seven  states  is  only  2,  as  com¬ 
pared  to  15  for  the  four  most  important  states. 


Thus  the  reactive  flux  is  not  simply  concentrated  into 
those  product  states  whose  internal  energies  are  closest  to 
the  internal  energy  of  the  reactant  state.  In  Table  111,  the 
internal  energy  and  the  state-selected  threshold  energy  are 
tabulated  for  initial  states  (0,11),  (0,12),  and  (0,13)  and, 
for  comparison,  for  selected  final  states  having  large  condi¬ 
tional  probabilities  for  reaction  out  of  these  three  initial 
states.  (The  four  largest  conditional  probabilities  for  each 
initial  state  are  listed  in  Table  1. )  For  those  final  states  corre¬ 
sponding  to  the  largest  conditional  probabilities,  there  is  a 
better  correlation  between  the  state-selected  thresholds,  de¬ 
fined  as  explained  above  in  terms  of  Eqs.  ( 4 )  and  ( 8 ) ,  for  the 
initial  and  final  states  than  between  the  internal  energies  of 
the  initial  and  final  states.  Note,  for  instance,  that  for  initial 
state  (0,11),  the  largest  conditional  probabilities  corre¬ 
sponding  to  final  states  in  the  v'  —  1  manifold  are  for  final 
states  (1,14)  and  (1,15)  (21%  and  16%,  respectively). 
These  are  the  two  final  states  whose  state-selected  thresholds 
( 1 .4 1 0  and  1 .484  e V ) ,  are  closest  to  the  state-selected  thresh¬ 
old  for  initial  state  ( 0, 1 1 )  (1 .484  eV ) .  The  internal  energies 
of  final  states  ( 1,14)  and  (1,15)  though,  1.293  and  1.306eV, 
are  considerably  higher  than  the  internal  energy  of  initial 
state  (0,11)  ( 1.157  eV).  The  two  final  states  in  the  v'  =  1 
vibrational  manifold  whose  internal  energies  are  closest  to 
the  internal  energy  of  initial  state  (0,11)  are  final  states 
(1,12)  and  (1,13).  However,  their  thresholds  are  much  low¬ 
er  in  energy  than  the  threshold  for  initial  state  ( 0, 1 1 )  and  the 
conditional  probabilities  Z}(1,12|0,11)  and  />(1,13|0,11) 
are  quite  small,  1%  and  3%. 

A  particularly  striking  example  of  this  trend  is  seen  for 
initial  state  (0,13).  The  largest  conditional  probability  cor¬ 
responding  to  a  final  state  in  the  v'  —  0  vibrational  manifold 
is  15%,  for  final  state  (0,21).  Initial  state  (0,13)  and  final 
state  (0,21)  have  nearly  identical  state-selected  thresholds 
(1.718  and  1.716  eV).  The  conditional  probability  corre¬ 
sponding  to  final  state  (0,22),  on  the  other  hand,  is  only  4%, 
even  though  initial  state  (0,13)  and  final  state  (0,22)  have 
nearly  the  same  internal  energies  ( 1 .454  and  1 .453  eV ) .  The 
threshold  for  final  state  ( 0,22 )  is  1 . 848  eV,  much  higher  than 
that  for  initial  state  (0,13). 

The  fact  that  the  threshold  energies  for  initial  and  final 
states  with  large  conditional  probabilities  are  in  better  agree¬ 
ment  than  the  internal  energies  of  the  states  suggests  that  the 
thresholds  are  correlated  with  transition  states  through 
which  the  reactive  flux  is  channeled.  This  view  is  in  agree¬ 
ment  with  the  work  in  Refs.  5  and  8,  where  it  is  shown  that 
quantized  transition  states  control  the  chemical  reactivity 
for  the  reaction  H  -f  Hj  —  -h  H. 

Tables  I  and  II  demonstrate  that  for  both  the  u  =  0  and  1 
vibrational  manifolds,  the  state-selected  threshold  reactivity 
is  primarily  vibrationaily  adiabatic  at  low  to  moderate  j 
(which  also  corresponds  to  low  to  moderate  energy);  i.e., 
the  important  final  states  tend  to  have  the  same  vibrational 
quantum  number  as  the  initial  state.  As  noted  above  for  se¬ 
lected  initial  states,  vibrational  nonadiabaticity  becomes 
more  important  at  higher  j  levels.  For  initial  states  with 
i’  —  0  and  j  =  1 1-14,  the  single  most  important  final  state 
has  u'  =  1;  for  initial  states  with  v  =  1  and  j  =  9-1 1,  the  sin¬ 
gle  most  important  final  state  has  v'  =  2. 
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TABLE  II.  Final-sute  specificity  of  reaction  out  of  initial  states  in  the  ti  =  I  manifold  for  the  O  -t-  H;  reaction  with  J  =  0. 


J 

P.K, 

(eV) 

(eV) 

25% 

50% 

75% 

inant  (v  ,f).  D(v‘ ,f\\.j)‘ 

N 

F 

N 

F 

N 

F 

Dom 

0 

0.17 

26 

2 

0.08 

3 

0.12 

5 

0.19 

(1,3)  19 

(1.4)  15 

(1.2)  16 

(1.5)  13 

1 

0.97 

0.17 

26 

2 

3 

0.12 

5 

0.19 

(1.3)  19 

(1.4)  19 

(1.2)  16 

(1.5)  14 

2 

0.97 

0.15 

26 

2 

0.08 

3 

0.12 

5 

0.19 

(1.3)  19 

(1,4)  19 

(1.2)  16 

(1.1)  14 

3 

0.98 

0.11 

26 

2 

0.08 

3 

0.12 

5 

0.19 

(1.4)  19 

(1.3)  19 

(1.2)  15 

(1.5)  14 

4 

1.14 

0.12 

32 

2 

0.06 

3 

0.09 

6 

0.19 

(1.8)  24 

(1.9)  19 

(1.7)  17 

(1.10)  7 

5 

1.22 

0.23 

37 

2 

0.05 

3 

0.08 

7 

0.19 

(1.10)  24 

(1. 11)  17 

(1.9)  17 

(1.8)  7 

6 

1.29 

0.22 

42 

2 

0.05 

3 

0.07 

7 

0.17 

(1.12)25 

(1.11)  22 

(1,7)  9 

(1.6)  6 

7 

1.38 

0.22 

48 

1 

0.02 

3 

0.06 

7 

0.15 

(1.13)  32 

(1.12)  13 

(1,14)  12 

(2,5)  7 

8 

1.47 

0.21 

52 

2 

0.04 

3 

0.06 

8 

0.15 

(1.14)  20 

(1.15)  15 

(2,8)  15 

(2.9)  11 

9 

1.56 

0.18 

56 

1 

0.02 

3 

0.05 

7 

0.13 

(2.10)  28 

(1.16)  17 

(1.15)  10 

(2.11)  9 

10 

1.66 

0.16 

65 

1 

0.02 

3 

0.05 

7 

0.11 

(2.12)  37 

(1. 17)  12 

(2,11)  9 

(2.9)  8 

11 

1.83 

0.21 

76 

1 

0.01 

3 

0.04 

9 

0.12 

(2.14)26 

(2.15)  16 

(1.19)  13 

(0,21)  5 

‘The  four  final  states  which  contribute  most  to  /‘f  and  their  normalized  conditional  probabilities  expressed  as  percentages,  X  100%,  listed  in 


order  of  the  size  of  D{v',f\l,j),  beginning  with  the  greatest  contributor. 


.V  Trends  in  vibrational  level  specificity  for  initial  states  in 
I  the  u  =  0  and  1  vibrational  manifolds  are  displayed  in  Figs.  3 
*  and  4.  For  initial  states  in  the  ground  vibrational  manifold, 
threshold  reactivity  is  vibrationally  adiabatic  up  to  J  =  5 
~  =  0.89  eV)  since  the  threshold  for  reaction  into  o  =  1 

is  0.97  eV.  Above  j  =  5,  the  fraction  of  reactivity  that  access¬ 
es  final  o'  =  1  states  increases  regularly,  until  by  y  =  1 1,  a 
greater  fraction  of  the  flux  goes  into  final  o'  =  1  sutes  than 
final  o'  =  0  states.  Reaction  into  o'  =  2  and  3  final  states 
remains  small  even  at  energies  considerably  above  the 
thresholds  for  reaction  into  these  vibrational  levels. 

A  similar  trend  is  seen  for  initial  states  in  the  o  =  1  vi¬ 
brational  manifold.  Threshold  reactivity  into  o'  =:  0  final 
states  remains  relatively  constant  over  the  entire  range  of  j 
values,  varying  between  8%  and  15%  of  the  reactive  flux. 
Beginning  with  y  =  6(1.29  eV),  reaction  into  final  o' =  2 


states  increases  regularly  until  j  =  10  and  then  drops  off  as 
reaction  begins  to  occur  into  o'  =  3  states.  By y  =  9,  a  greater 
fraction  of  the  flux  out  of  initial  o  =  1  states  goes  into  final 
o'  =  2  states  than  into  o'  =  1  states. 

It  is  seen  in  Tables  I  and  II  that  initial  states  within  a 
given  vibrational  manifold  whose  thresholds  occur  at  nearly 
the  same  energy  [(0,0),  (0,1),  (0,2),  (0,3),  and  (0,4); 
(0,5),  (0,6),  and  (0,7);  or  (1,0),  (1,1),  (1,2),  and  (1,3)] 
send  flux  to  the  same  or  nearly  the  same  set  of  final  states. 
This  phenomenon  may  be  due  to  the  dominance  of  the  indi¬ 
vidual  threshold  reactivities  of  a  set  of  initial  states  by  the 
same  energy  level  of  the  transition  state.®  Like  other  reson¬ 
ances,  quantized  transition  state  thresholds  are  interme¬ 
diates,  and  so  the  threshold  flux  through  a  quantized  transi¬ 
tion  state  tends  to  lose  memory  of  its  particular  state  of 
origin. 
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Hi  B.  Channel-selected  specificities  for  H-i-H} 

To  learn  whether  these  trends,  which  we  discussed 
above  for  O  -I-  H2 ,  are  an  isolated  phenomenon  for  that  case 
or  whether  they  might  be  more  general,  we  present  a  similar 


TABLE  III.  Internal  energies  E,„,  and  state-selected  threshold  energies 
in  electron-volts  for  selected  initial  and  final  states  of  the  O  -f  H,  reaction. 


Initial  states 

Final  states 

V 

J 

^tni  ^ihr 

y' 

/ 

■^ini  ^ihr 

v’ 

/ 

^thr 

0 

11 

1.157  1.484 

0 

17 

1.031  1.358 

1 

12 

1.129  1.272 

0 

12 

1.302  1.596 

0 

18 

1.108  1.396 

1 

13 

1.184  1.350 

0 

13 

1.454  1.718 

0 

19 

1.189  1.510 

1 

14 

1.243  1.410 

0 

20 

1.274  1.594 

1 

15 

1.306  1.484 

0 

21 

1.362  1.716 

1 

16 

1.373  1.568 

0 

22 

1.453  1.848 

1 

17 

1.444  1.646 

1 

18 

1.516  1.730 

1 

19 

1.596  1.844 

FIG.  2.  Normalized  product  distribution  P(i;',/|l,10)  and  state-to-state 
reaction  probabilities  vs  final  rotational  quantum  number  f  for  the 
0-k  Hj  reaction  withy  =0  for  initial  state  (1,10)  at  energy  P.h,  =  1.66eV. 
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FIG.  3.  Normalized  product  distributions  D(v'\0,j),  which  indicate  the 
percentage  of  threshold  reactive  flux  into  each  final  vibrational  level  ti'  vs 
initial  rotational  quantum  number y  for  initial  states  in  the  ground  vibration¬ 
al  manifold  at  their  threshold  energies  for  the  O  -l-  reaction  with  J  =  0. 


FIG.  4.  Normalized  product  distributions  D{v'\\,j),  which  indicate  the 
percentage  of  threshold  reactive  flux  into  each  final  vibrational  level  u'  vs 
initial  rotational  quantum  number  j  for  initial  states  with  ti  =  I  at  their 
threshold  energies  for  the  O  Hj  reaction  with  7  =  0. 


Study  for  the  H  -f  Hj  reaction.  Results  at  the  fully  specified 
channel  level  are  summarized  in  Tables  IV-VI.  For  J  =  0, 
where  channels  and  states  are  the  same,  three  initial  states 
with  different  reaction  thresholds  were  arbitrarily  selected 
for  study  { (0,4),  ( 1,4),  and  (0,10)  ].  It  is  seen  in  Table  III 
that  threshold  flux  for  these  states  is  very  final-state  specific; 
in  each  case,  a  single  final  state  accounts  for  over  23%  of  the 
reactive  flux,  two  states  account  for  over  50%,  and  at  most 
five  are  needed  to  account  for  75%  of  the  flux.  The  largest 
state-to-state  reaction  probabilities  are  all  vibrationally  adia¬ 
batic. 

In  order  to  test  whether  these  trends  also  prevail  for 
nonzero  values  of  the  total  angular  momentum,  we  evaluat¬ 
ed  results  of  the  H  -f  Hj  reaction  at  /  =  4.  Since  for  nonzero 
J  both  values  of  the  parity  are  possible,  results  for  /  =  4, 
F  =  -f-  1  and  for  7  =  4,  F  =  —  1  are  summarized  separately 
in  Tables  V  and  VI,  respectively. 

There  are  many  more  channels  for  7  =  4  than  for  7  =  0, 
since  /  can  vary  from  min(/  —  4,0)  toy  4  for  a  given  value 
of  j.  Nevertheless,  the  value  of  N(.x% )  remains  small;  the 
flux  out  of  initial  channels  is  quite  final-channel  specific.  For 
example.  Table  VII  shows  channel-to-channel  reaction 
probabilities  for  initial  channels  (0,10,6)  and  (0,10,8)  and 


final  channels  (0,8,/ ')  and  (0,9,/ ')  for  all  allowed  values  of 
the  final-channel  orbital  angular  momentum  quantum  num¬ 
ber  / '  for  7F  =  4  -f  .  For  a  given  initial  channel  and  final 
state,  one  of  the  channel-to-channel  reaction  probabilities  is 
larger  than  the  others  by  more  than  an  order  of  magnitude. 
For  example,  for  initial  channel  (0,10,6),  the  channel-to- 
channel  reaction  probability  into  final  channel  (0,8,12), 
^fo.io.6)(o.i.i2).is0106,butFfo,o,6„o.8.;)  is  less  than  0.01  for 
/ '  =  4, 6, 8,  or  10.  Similarly,  the  largest  value  of  / '  ( 1 3 )  domi¬ 
nates  reaction  out  of  initial  channel  (0,10,6)  into  final  state 
(0,9 ) .  The  data  for  other  final  states  also  show  that  reactants 
in  initial  channel  (0,10,6)  react  preferentially  to  form  prod¬ 
ucts  with  the  highest  allowed  value  of  / '  for  a  given  v'  and y'. 
Table  VII  shows  that  for  initial  channel  (0,10,8),  the  second 
highest  value  of  / '  for  a  given  v'  and y'  is  preferred.  A  similar 
trend  is  found  for  the  other  initial  channels  (0,10,/)  with 
different  values  of  /,  so  that  a  general  summary  of  this 
propensity  for  (0,10,/)  initial  channels  would  be 
y,  —  /,  =  —  ( yy  —  //).  where  i  and /  denote  initial  and  final. 

Hi  C.  State-selected  specificities  for  H+H, 

As  mentioned  above,  state-selected  and  channel-select¬ 
ed  specificities  are  the  same  for  7  =  0.  Although  our  primary 


TABLE  IV.  Final-state  specificity  of  reaction  out  of  selected  initial  states  for  the  H  H,  reaction  with  7=  0. 


(eV) 

^rstibr 

(cV) 

25% 

50% 

75% 

\o,jr 

N 

F 

N 

F 

N 

F 

Dominant  (u’,/).  D(v',f\ 

(0,4) 

14 

1 

0.07 

2 

0.14 

4 

0.29 

(0,4)  31  (0,3)  21  (0,5)  17 

(0,0)  12 

(1.4) 

1.18 

0.27 

20 

1 

0.05 

2 

5 

0.25 

(1,3)  33  (1,4)  24  (1,1)  6 

(0,0)  5 

(0,10) 

1.49 

0.48 

31 

I 

0.03 

2 

0.06 

4 

0.13 

(0,10)  37  (0,9)  29  (0,8)  6 

(1.7)  5 

*The  four  final  states  which  contribute  most  to  F*  and  their  normalized  conditional  probabilities  expressed  as  percenuges,  /7(v',/|i;,;)  X  100%,  listed  in 
order  of  the  size  of  f>(u',/|v,y),  beginning  with  the  greatest  contributor. 
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TABLE  V.  Final-channel  specificity  of  reaction  out  of  selected  initial  channels  for  the  H  -I-  Hj  reaction  with  JP=i  +  . 


IvJ.l) 

(eV) 

E 

^rtUhr 

(eV) 

25% 

50% 

75% 

75% 

Dominant  iv'.fj'). 

N„p„ 

N 

F 

V 

F 

N 

F 

(0.4,0) 

0.900 

0.50 

51 

2 

0.04 

4 

0.08 

8 

0.16 

(0,2.6)  19 

(0,1,5)  16 

(0,3.7)  14 

(0,0,4) 

8 

(0.4.2) 

0.900 

0.50 

51 

2 

0.04 

3 

0.06 

7 

0.14 

(0.2.6)  19 

(0,1.5)  18 

(0.3.7)  13 

(0,0.4) 

11 

(0,4,4) 

0.900 

0.50 

51 

2 

0.04 

5 

0.10 

9 

0.18 

(0,3.5)  15 

(0,4.6)  14 

(0.4,4)  12 

(0.5,5) 

8 

(0,4,6) 

0.900 

0.50 

51 

2 

0.04 

4 

0.08 

10 

0.20 

(0.4,4)  15 

(0.3,5)  13 

(0.4.6)  12 

(0.5.5) 

10 

(0.4.8) 

o.roo 

0.50 

51 

2 

0.04 

5 

0.10 

10 

0.20 

(0,5,1)  15 

(0,4.0)  12 

(0.3,1)  9 

(0.6.2) 

8 

(1.4.0) 

1.205 

0.29 

81 

3 

0.04 

5 

0.06 

11 

0.14 

(1.2.6)  12 

(1,1,5)  12 

(1.3,1)  12 

( 1.4.0) 

8 

(1.4.2) 

1.205 

0.29 

81 

2 

0.03 

5 

0.06 

13 

0.16 

(1.1.5)  13 

(1.2.6)  13 

(1,0,4)  9 

(1.3.7) 

8 

(1.4.4) 

1.205 

0.29 

81 

2 

0.03 

5 

0.06 

13 

0.16 

(1.3.5)  16 

(1,4,4)  10 

(1.3.3)  9 

(1.4,6) 

9 

(1.4.6) 

1.205 

0.29 

81 

2 

0.03 

5 

0.06 

IS 

0.19 

(1,3.5)  16 

(1.4.4)  II 

(1.3,3)  9 

(1,2,4) 

8 

(1.4,8) 

1.205 

0.29 

81 

3 

0.04 

7 

0.09 

14 

?.17 

(1.2,6)  11 

(0,8.4)  9 

00 

(1.1.5) 

7 

(0.10,6) 

1.500 

0.49 

126 

1 

0.01 

3 

0.02 

5 

0.04 

(0,8,12)  31 

(0,9.13)  19 

(0,7,11)  17 

(1.6.10) 

7 

(0,10.8) 

1.500 

0.49 

126 

1 

0.01 

2 

0.02 

6 

0.05 

(0,9,11)  29 

(0,8,10)  23 

(0,10,12)  9 

(1,6,8) 

5 

(0,10,10) 

1.300 

0.49 

126 

2 

0.02 

3 

0.02 

9 

0.07 

(0,9.9)  25 

(0,10,10)  23 

(0.8,8)  10 

(0,9,11) 

5 

(0,10.12) 

1.500 

0.49 

126 

1 

0.01 

2 

0.02 

8 

0.06 

(0,10,8)  32 

(0.9.7)  18 

(0.11,9)  8 

(0,10.10) 

4 

(0.10,14) 

1.500 

0.49 

126 

1 

0.01 

2 

0.02 

6 

0.05 

(0,10.6)  30 

(0,11.7)  27 

(0.9,5)  8 

(1.8.4) 

4 

‘The  four  final  channels  which  contribute  most  to  /*  f  and  their  normalized  conditional  probabilities  expressed  as  percentages,  />(  v',/,/ '  |  x  100%,  listed 

in  order  of  the  size  of  D(  u',/,l  beginning  with  the  greatest  contributor. 


TABLE  VI.  Final-channel  specificity  of  reaction  out  of  selected  initial  channels  for  the  H  -l-  H,  reaction  with  JP  =  A  —  . 


(U,/,/) 

(eV) 

^rel.thr 

(cV) 

25% 

50% 

75% 

Dominant  (v',/,/'),  D(v',/J'\v,J,ly 

iV 

F 

S 

F 

N 

F 

(0.4.1) 

0.900 

0.50 

37 

2 

0.05 

5 

0.14 

8 

0.22 

(0.3.6)  16 

(0,2,5)  16 

(0,4,7)  9 

(0,3.6)  8 

(0.4,3) 

0.900 

0.50 

37 

2 

0.05 

4 

0.11 

7 

0.19 

(0.2,5)  17 

(0,3,6)  15 

(0,2,3)  11 

(0.1.4)  10 

(0.4.5) 

0.900 

0.50 

37 

2 

0.05 

4 

0.11 

6 

0.16 

(0.2,5)  17 

(0,2,3)  17 

(0,1,4)  14 

(0.3.4)  11 

(0,4.7) 

0.900 

0.50 

37 

3 

0.08 

5 

0.14 

9 

0.24 

(0,2,5)  11 

(0,4,3)  11 

(0.3,2)  11 

(0.2.3)  9 

(1.4,1) 

1.205 

0.29 

61 

2 

0.03 

5 

0.08 

9 

0.15 

(1.3,2)  15 

(1.2,3)  14 

(1.2.5)  9 

(1.3.6)  9 

(1.4.3) 

1.205 

0.29 

61 

2 

0.03 

4 

0.07 

8 

0.13 

(1.2,3)  18 

(1.3.2)  15 

(1,1.4)  10 

(1.2.5)  9 

(1,4.5) 

1.205 

0.29 

61 

2 

0.03 

3 

0.05 

9 

0.15 

(1,2,3)  18 

(1.1.4)  18 

(1.2,5)  15 

(1.3.2)  8 

(1.4.7) 

1.205 

0.29 

61 

2 

0.03 

5 

0.08 

12 

0.20 

(1.2,5)  17 

(1.1.4)  11 

(1.3,6)  9 

(0.8.5)  7 

(0,10,7) 

1.500 

0.49 

95 

1 

0.01 

2 

0.02 

5 

0.05 

(0,8,1 1)  29 

(0,9,12)  25 

(0.7,10)  10 

(1.6,9)  7 

(0,10.9) 

1.500 

0.49 

95 

1 

0.01 

3 

0.03 

8 

0.08 

(0,9,10)  29 

(0,8,9)  16 

(0,10,11)  14 

vO 

00 

1“^ 

(0.10.11) 

1.500 

0.49 

95 

1 

0.01 

3 

0.03 

10 

0.11 

(0,10,9)  26 

(0,9,8)  22 

(0,9.10)  5 

(0,10.11)  5 

(0.10,13) 

1.500 

0.49 

95 

1 

0.01 

3 

0.03 

8 

0.08 

(0,10,7)  33 

(0,11,8)  16 

(0,9.6)  13 

(1.8,5)  4 

‘The  four  final  channels  which  contribute  most  to/’*  and  their  normalized  conditional  probabilities  expressed  as  percentages.  D(.  u\j,l’\u,j,l)  x  100%,  listed 
in  order  of  the  size  of  />(u',/,/'|i;,y,/),  beginning  with  the  greatest  contributor. 


TABLE  VII.  Channel-to-channel  reaction  probabilities  for  initial  channels 
(0,10,6)  and  (0,10,8)  and  final  channels  (0,8,/')  and  (0,9,/')  for  all  allowed 
values  of  the  final-channel  orbital  angular  momentum  quantum  number  /' 
for  the  H  -h  H2  reaction  with  JP=i  +  . 


1/ 

/ 

/' 

pM  • 

* 

pK  » 

•  (aio.i)it 

0 

8 

4 

5.28(  -4)'’ 

7.62(  -  6) 

0 

8 

6 

4.29(  -  6) 

4.44(  -  4) 

0 

8 

8 

1.30(  -  5) 

2.28(  -  3) 

0 

8 

10 

1.32(-3) 

5.41(-2) 

0 

8 

12 

1.06(  -  1) 

2.44(  -  3) 

0 

9 

5 

8.96(  -  5) 

1.07(  -6) 

0 

9 

7 

3.06(  -  6) 

9.35(  -  5) 

0 

9 

9 

8.03(  -  6)  ' 

1.55(  -3) 

0 

9 

11 

8.22(  -4) 

6.79(  -  2) 

0 

9 

13 

6.31(  -2) 

2.16(  -3) 

indicates  channel-to-channel  reaction  probability  F*,.  for  initial 
channel  n  = 

"  Numbers  in  parentheses  are  powers  of  ten. 


emphasis  here  is  on  the  completely  resolved  channel-select¬ 
ed  level  of  detail,  Table  VIII  presents  the  /  =  4  results  for 
H  -F  Hj  in  a  state-selected  way.  Here  up  to  nine  /  values 
contribute  to  the  results  for  a  given  initial  state,  and  up  to 
nine  final  / '  values  are  included  for  a  given  final  state. 

For  /  =  4,  as  for  /  =  0,  threshold  reactivity,  even  when 
averaged  over  /  values,  is  highly  final-state  specific.  At  most 
two  final  states  are  required  to  account  for  25%  of  the  flux, 
and  about  five  states  account  for  over  75%  of  the  flux.  To  a 
large  extent,  the  same  final  states  dominate  state-selected 
threshold  reactivity  for  both  values  of  J.  For  initial  states 
(0,4)  and  (1,4),  three  of  the  four  final  states  that  contribute 
most  to  P  ^  are  the  same  as  for  /  =  0,  although  their  order¬ 
ing  according  to  relative  importance  may  differ.  Two  of  the 
four  dominant  final  states  for  initial  state  (0,10)  are  com¬ 
mon  to  both  J  =  0  and  4. 

The  7  =  0  reaction  tends  to  be  somewhat  more  final- 
state  specific  than  the  7  =  4  reaction,  as  evidenced  by  the 


4320 


Chattietd,  Truhlar,  and  Schwenke:  State-selected  reaction  dynamics 


TABLE  VIIl.  Final-state  specificity  of  reaction  out  of  selected  initial  states  for  the  H  -i-  H;  reaction  with  y  =  4. 


{u,j) 

(eV) 

(eV) 

23% 

50% 

73% 

Dominant  (v'.j).  D(v.j  \i\j)‘ 

N 

F 

N 

F 

N 

F 

(0,4) 

0.900 

0.30 

14 

2 

0.14 

3 

0.29 

(0.3)  25 

(0.2)  21 

(0,4)  21 

(0,1)  13 

(1,4) 

1.203 

0.29 

KB 

2 

3 

0.30 

(1.3)  23 

(1.2)  19 

(1,4)  15 

(1,1)  11 

(0,10) 

1.300 

0.49 

D 

1 

3 

0.16 

(0.9)  25 

(0.8)  23 

(0,10)  14 

(0,7)  8 

‘The  four  final  states  which  contribute  most  to  /•*  and  their  normalized  conditional  probabilities  expressed  as  percentages,  D{v'.j  \u,j)  x  100%,  listed  in 
order  of  the  size  of  D{v'./\u,j),  beginning  with  the  greatest  contnbutor. 


slightly  lower  values  of  N{x% )  and  slightly  higher  values  of 
Div',f\vJ)  for  /=  0  than  for  7=  4.  The  vibrational  level 
specificity  of  the  7  =  0  and  4  reactions  is  the  same  for  initial 
states  (0,4)  and  ( 1,4)  [  100%  of  the  flux  is  into  u'  =  0  for 
(0,4);  22%  is  into  u'  =  0  and  78%  is  into  v'  =  1  for  ( 1,4)  J. 
Vibrationally  nonadiabatic  processes  account  for  a  slightly 
larger  fraction  of  the  reactive  flux  for  initial  state  (0,10) 
when  7=4  than  when  7  =  0  ( 86%,  10%,  and  4%  of  the  flux 
is  into  v'  =  0,  1,  and  2,  respectively,  for  7=  0;  82%,  18%, 
and  0%  is  into  v'  =  0,  1,  and  2  for  7  =  4). 

Ill  O.  Above  threshold  specificity 

Although  the  main  emphasis  of  the  present  study  is  on 
specificity  at  threshold,  we  also  investigated  briefly,  for  com¬ 
parison,  how  the  final-state  specificity  of  state-selected  reac¬ 
tivity  changes  as  energy  is  increased  above  threshold.  Table 
IX  summarizes  such  a  study  for  initial  state  (0,5)  of  the 
O  -I-  Hj  reaction  with  7  =  0  in  the  energy  range  0.70-1.90 
eV.  As  the  energy  is  increased  above  threshold  (0.89  eV), 
the  state-selected  reactivity  becomes  less  final-state  specific. 
By  1.90  eV,  the  four  most  important  final  states  account  for 
less  than  25%  of  the  total  reactive  flux,  and  the  22  most 
important  final  states  account  for  less  than  75%  of  the  flux, 
whereas  at  the  threshold  energy,  two  final  states  contribute 
more  than  25%  of  the  flux,  and  four  final  states  contribute 


more  than  75%  of  the  flux.  This  behavior  of  initial  state 
(0,5)  at  high  energy  also  stands  in  contrast  to  threshold 
reactivity  summarized  in  Tables  I  and  II.  For  example,  at  the 
highest-energy  threshold  studied,  that  for  initial  state  (0,14) 
at  1.86  eV,  over  25%  of  the  state-selected  flux  goes  into  a 
single  final  state  (1,19). 

In  Figs.  5(a)  and  5(b),  N(x%)  and  F(x%)  are  plotted 
vs  the  total  energy  at  threshold.  Figure  5(a)  indicates  that 
final-state  specificity  is  greatest  at  the  state-selected  reaction 
threshold.  Both  above  and  below  threshold,  the  reaction  is 
less  specific.  [Even  at  energies  below  the  threshold  for  initial 
state  (0,5),  there  is  still  a  small  probability  of  reaction  ] 
Above  about  0.90  eV,  final-state  specificity  as  indicated  by 
N(x%)  and  D(v',f  \0,5)  (Table  IX)  decreases  markedly; 
the  number  of  small  contributors  increases  dramatically 
with  energy,  as  indicated  by  the  large  value  of  N(15%)  at 
high  energies.  Although  the  specificity  docs  decrease  at  ener¬ 
gies  above  threshold,  some  specificity  remains.  For  example, 
five  final  channels  still  account  for  25%  of  the  flux  at  1.90 
eV,  and  Fig.  5(b)  emphasizes  that  the  fraction  of  final  states 
(F(x%)]  that  account  for  a  given  percentage  of  the  flux 
does  not  necessarily  increase  with  energy  above  threshold. 
At  1.90  eV,  the  fraction  is  0.06  for  25%  of  the  flux,  as  com¬ 
pared  to  0.08  at  threshold. 

Up  to  1.90  eV,  the  four  largest  state-to-state  reaction 
probabilities  for  reaction  out  of  (0,5)  are  all  vibrationally 


TABLE  IX.  Final-state  specificity  of  reaction  out  of  initial  state  (u  =  0,j  =  S)  for  the  O  -i-  Hj  reaction  with  7  =  0. 


(eV) 

(eV) 

25% 

50% 

75% 

ninant  (v,f),  Z)(u',/|0,5)' 

F 

S 

F 

S 

F 

Dor 

0.70 

MM 

12 

2 

0.17 

4 

0.33 

6 

0.30 

(0.6)  17 

(0,5)  16 

(0,7)  15 

(0,4)  13 

16 

2 

0.13 

3 

0.19 

3 

0.31 

(0,8)  22 

(0,9)  22 

(0.7)  16 

(0,10)  14 

0.90 

0.41 

23 

2 

0.09 

3 

0.13 

4 

0.17 

(0,10)  23 

(0.9)  23 

(0,11)  16 

(0,8)  15 

1.00 

0.31 

27 

2 

0.07 

3 

0.11 

6 

0.22 

(0,10)  20 

(0.1 1)  19 

(0,9)  14 

(0.12)  12 

1.10 

0.61 

30 

2 

0.07 

4 

0.13 

9 

0.30 

(0.11)  18 

(0,10)  16 

(0,12)  12 

(0.9)  10 

1.20 

0.71 

34 

2 

0.06 

5 

0,15 

II 

0.32 

(0,11)  16 

(0.12)  13 

(0,10)  12 

(0.9)  5 

1.30 

0.81 

43 

3 

0.07 

8 

0.19 

15 

0.35 

(0.11)  11 

(0.12)  9 

(0,10)  7 

(0,5)  5 

1.40 

0.91 

49 

4 

0.08 

10 

0.20 

17 

0.33 

(0.11)  7 

(0,15)  7 

(0,12)  6 

(0,14)  6 

1.30 

1.01 

33 

4 

0.08 

10 

0.19 

17 

0.32 

(0.15)  8 

(0,14)  8 

(0,2)  6 

(0.3)  5 

1.60 

l.II 

38 

3 

0.05 

8 

0.14 

18 

0,31 

(0.15)  9 

(0,14)  8 

(0.2)  8 

(0.1)  7 

1.70 

1.21 

68 

3 

0.04 

8 

0.12 

18 

0.26 

(0,13)  9 

(0,2)  9 

(0.14)  8 

(0.1)  8 

1.80 

1.31 

76 

4 

0.05 

10 

0.13 

21 

0.28 

(0.2)  8 

(0.1)  8 

(0,15)  6 

(0,14)  6 

1.90 

1.41 

79 

3 

0.06 

11 

0.14 

23 

0.29 

(0,2)  8 

(0,1)  8 

(0.17)  6 

(0,15)  6 

‘The  four  final  states  which  contribute  most  loP^  and  their  normalized  conditional  probabilities  expressed  as  percentages,  D{  i;',/|0,S )  x  IOO%,  listed  in 
order  of  the  size  of  />(u',/|0,3),  beginning  with  the  greatest  contributor. 
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Energy  (eV) 


Energy  (eV) 

FIG.  5.  For  the  O  -(-  Hj  reaction  with  7  =  0  (a)  the  minimum  number 
S(x% )  of  sute-to-sute  reaction  probabilities  needed  to  account  for  x%  of 
the  flux  out  of  initial  state  (0,3)  vs  energy  for  x  =  23,  30,  and  73;  (b)  the 
fraction  F(x%)  =  N(x%)/N^  vs  energy  for  x  =  23,  30,  and  73. 


adiabatic,  in  contrast  to  the  high  degree  of  vibrational  nona- 
diabaticity  for  the  (0,14)  and  other  high-y  states  at  similar 
energies.  At  1.86  eV,  57%  of  the  reactive  transitions  out  of 
(0,5)  are  vibrationally  adiabatic,  compared  to  only  21%  out 
of  (0,14). 

Figure  6  shows  how  the  product  distribution  specifici¬ 
ties  D{v',f\v,j)  change  with  energy  for  initial  state  (0,5) 
and  v'  =  0,  which  is  the  dominant  final  vibrational  level.  At 
0.7  eV,  which  is  below  the  state-selected  threshold  of  0.89 
eV,  i)(0,/|0,5)  as  a  function  of f  is  a  single,  relatively  broad 
peak  with  a  maximum  at /  =  6.  By  0.9  eV,  very  close  to  the 
threshold  energy,  this  peak  has  shifted  to  the  right  (the  max¬ 
imum  is  now  at/  =10)  and  grown  both  taller  and  narrower. 
At  the  same  time,  a  second  peak  has.begun  to  appear  at  low 
/.  This  bimodal  distribution  is  maintained  as  the  energy  is 
increased  further.  The  peak  at  low  f  becomes  more  pro¬ 
nounced  with  increasing  energy,  while  the  peak  at  high  f 
shifts  to  the  right  and  grows  broader.  By  1.9  eV,  the  highest 
peak  is  the  one  at  low  f.  These  trends  are  reflected  in  the 


0  5  10  15  20  25 

j' 

FIG.  6.  Normalized  product  distributions  /)(i)',/|0.3)  vs  final  rotational 
quantum  number  /  for  initial  stale  (0.3)  and  Anal  vibrational  quantum 
number  u'  equal  to  0  for  five  energies  I'or  the  O  -I-  H,  reaction  with  7  =  0. 
The  energies  range  from  0.2  eV  below  the  state-selected  threshold  to  1 .0  eV 
above  it. 


dominant  final  states  and  their  conditional  probabilities  giv¬ 
en  in  Table  IX.  The  broadening  of  the  high-/  peak  above 
threshold  is  reflected  in  the  decreasing  conditional  probabili¬ 
ties,  and  the  increasing  bimodality  of  the  distribution  is  evi¬ 
denced  by  the  fact  that  by  1.8  eV,  (0,2)  and  (0,1)  arc  the 
most  important  final  states.  Since  the  highest  and  most  con¬ 
centrated  peak  occurs  at  0.9  eV,  Fig.  6  supports  the  observa¬ 
tion  made  above  that  state-selected  reactivity  is  most  specific 
at  a  state-selected  threshold. 

Another  way  to  study  trends  in  final-state  specificity  as 
energy  is  increased  beyond  threshold  is  to  choose  one  energy 
and  evaluate  the  final-state  specificity  of  several  states  whose 
thresholds  occur  over  a  range  of  energies  up  to  the  chosen 
energy.  This  strategy  has  the  advantage  that  the  number  of 
open  final  states  remains  constant.  The  results  of  such  a 
study  are  depicted  in  Fig.  7  for  initial  states  in  the  ground 
vibrational  manifold  at  1.86  eV.  For  initial  states  not  too  far 
above  threshold  (y  =  9-14),  A^(jc%)  decreases  regularly  as 
the  threshold  is  approached  until  j  =  14,  for  which  1 .86  eV  is 
the  threshold  energy.  For  j=  15,  whose  threshold  is  above 
1.86  eV,  the  reaction  is  less  specific  than  for  j  =  14.  Thus  we 
find  that  for  a  given  energy  E,  the  final-state  specificity  is 
greatest  for  those  initial  states  whose  thresholds  are  at  or 
very  close  to  E,  just  as  before  we  found  that  for  a  given  initial 
state,  the  specificity  is  greatest  at  its  threshold  energy.  This 
again  confirms  our  conclusion  that  final-state  specificity  is 
especially  large  at  state-selected  reaction  thresholds. 

In  Fig.  8,  the  fraction  of  reactive  flux  into  each  final 
vibrational  level  at  1 . 86  e  V,  given  by  the  product  distribution 
D(v'\0,j),  is  plotted  as  a  function  of  j  for  initial  states  in  the 
ground  vibrational  manifold.  For  most  of  the  range  of  j  val¬ 
ues,  in  particular  for  j  =  3-12,  we  find  that  50%-60%  of  the 
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FIG.  7.  The  minimuin  number  N(x% )  of  state-to-state  reaction  probabili¬ 
ties  needed  to  account  for  x%  of  the  reactive  flux  and  fraction 
F(x%'i  =  N(x%)/N^^  vs  initial  rotational  quantum  numbery  for  lilial 
-'ites  in  the  ground  vibrational  manifold  for  x  =  25,  50,  and  75  for  the 
O  -i-  H,  reaction  with  total  energy  1.86  eV. 


reactive  flux  is  vibrationally  adiabatic.  At  low  values  of  j 
{ 7  =  0-2 ) ,  the  degree  of  adiabaticity  is  more  en  atic,  and  for 
high  j  (j  =  13  and  14),  whose  thresholds  are  near  1.86  eV, 
the  degree  of  nonadiabaticity  increases  substantially. 

IV.  CONCLUDING  REMARKS 

A  high  degree  of  final-state  specificity  has  been  demon¬ 
strated  for  sute-selected  reactivity  at  energies  near  the  state- 
selected  reaction  thresholds  for  the  reactions 
O  -t-  —OH  +  H  and  H  -I-  Hj  —  -f  H.  This  specificity 
has  been  found  for  nonzero  values  of  the  total  angular  mo¬ 
mentum  J  as  well  as  for  y  =  0,  and  it  is  particularly  marked 
for  large  values  of  the  rotational  quantum  number. 

We  do  not  know  of  any  previous  study  that  has  demon- 


F1G.  8.  Normalized  product  distributions  D{v'\Q.j),  which  indicate  the 
percentage  of  threshold  reactive  flux  into  each  final  vibrational  level  v  vs 
initial  rotational  quantum  numbery  for  initial  states  in  the  ground  vibration¬ 
al  manifoldata  total  energy  of  1.86  eV  for  the  O  -t-  H,  reaction  withy  =  0. 


strated  the  remarkable  specifici'ies  that  emerge  from  these 
accurate  quantum  dynamics  .alculations.  The  present  re¬ 
sults  are  particularly  dramatic  for  O  +  H,  because  the  final 
diatomic,  OH,  is  not  the  same  as  the  initial  one.  H, ,  and  yet 
the  specificity  is  very  strong.  One  might  have  been  less  sur- 
pnsed  if  such  specificities  were  observed  only  in  symmetnc 
reactions  such  as  H  4-  H, .  The  closest  previously  observed 
analog  of  this  specificity  seems  to  be  associated  with  inelastic 
processes  of  the  form +BC(u,j)—A  -t- flCtp'.y'', ,  where 
indeed  the  initial  and  final  diatomic  are  the  same.  For  exam¬ 
ple,  trajectory  studies  of  several  atom-diatom  systems  have 
shown  a  strong  tendency  for  tutational  energy  changes  to 
compensate  viSrational  energy  changes.'^  ”  High  rotational 
levels  were  found  to  enhance  vibrationally  inelastic  cross 
sections,’’'"’  and  the  results  show  that  the  degree  of  final- 
state  specificity  increases  with  initial  J  for  both  vibrationally 
elastic  and  inelastic  processes."  '*  ”  The  most  dramatic  ex¬ 
perimental  examples  are  found  in  the  work  of  Pntchard  and 
CO- workers”*®  for  rare  gas  collisions  with  Lij.  There  the 
effect  Mas  seen  for  high  initial  j  and  low  relative  translational 
energies  and  it  manifested  itself,  as  in  the  trajectory  studies, 
as  a  narrow  distribution  of  final  j  values  for  each  final  v'. 
Although  those  collisions  are  quite  different  from  the  reac¬ 
tive  ones  studied  here,  they  might  provide  useful  analogies  to 
pursue  in  future  work. 

Another  possible  avenue  of  future  study  is  to  relate  the 
final-state  distributions  to  the  properties  of  the  transition 
state.  If  the  transition  state  is  viewed  as  a  short-lived  inter¬ 
mediate,***'  this  viewpoint  is  similar  to  the  transition- 
state-product-state  overlap  model  of  direct  bimolecular  re¬ 
actions,**  or  to  the  interpretation  of  product  rotational  state 
distributions  in  photodissociation  processes  in  terms  of  the 
properties  of  initial  dissociative  resonance  states.**  We  have 
already  applied  this  kind  of  analysis  to  some  aspects  of  the 
product  rotational  distributions  for  the  H  -|-  Hj  reaction,* 
and  we  believe  that  this  kind  of  analysis  should  also  be  useful 
for  understanding  .he  high  state-to-state  specificity  of  near¬ 
threshold  processes  observed  in  the  present  study. 
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